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摘 ”要 : 合成 孔径 雷达 (Synthetic Aperture Radar, 简称 SAR) 图 像 受 到 斑点 噪声 影响 以 及 成 像 条 件 的 
变化 ， 使 得 同一 场景 的 两 幅 SAR 图 像 之 间 存 在 很 大 差异 ， 利 用 基于 边缘 特征 或 灰 度 信息 进 
行 SAR 图 像 配 准 ， 难 以 达到 预期 效果 . 基于 此 ， 本 文 提出 一 种 基于 分 割 区 域 的 SAR 图 像 配 准 方 
法 ， 该 方法 首先 利用 主动 轮 廊 方 法 对 去 噪 后 的 SAR 图 像 进行 分 割 ， 得 到 图 像 分 割 区 域 ， 然 后 提 
取 分 割 区 域 的 特征 ， 利 用 区 域 匹 配 度 函 数 对 图 像 进行 粗 匹 配 ; 最 后 利用 互信 息 优 化 粗 匹 配 下 的 变 
换 参 数 ， 实 现 由 粗 到 细 的 图 像 配 准 . 实验 结果 验证 了 算法 的 有 效 性 . 
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图 像 配 准 是 将 不 同时 间 、 不 同 角 度 或 不 同 传感器 所 得 的 两 幅 或 两 幅 以 上 的 图 像 进行 几何 对 
准 的 处 理 过 程 43， 方 法 主要 有 基于 区 域 和 基于 特征 的 图 像 配 准 . 

基于 区 域 的 图 像 配 准 方法 有 空间 相关 法 加、 不 变 矩 法 和 频 域 相关 法 向 等， 这些 方法 直接 利 
用 图 像 的 灰 度 信息 ， 配 准 精 度 高 ， 但 是 对 图 像 灰 度 变化 比较 敏感 ， 计 算 复 杂 度 高 ， 对 目标 的 旋 
转 等 比较 敏感 ， 因 而 这 些 方法 不 能 适应 SAR 图 像 配 准 . 近 几 年 ， 基 于 互信 息 最 大 化 的 相似 度 
度量 在 SAR 图 像 配 准 中 得 到 了 应 用 .Xieig 直接 利用 互信 息 方 法 将 去 噪 后 的 SAR 图 像 配 准 ， 
去 品 效 果 影 响 算法 的 鲁 棒 性 且 计 算 量 大 .Hiepl6l 提出 了 一 种 自动 分 层 的 SAR 图 像 配 准 方法 ， 
该 方法 无 需 对 SAR 图 像 进行 去 品 ， 但 采用 插值 作为 搜索 策略 计算 杂 度 大 .Lil7',8l 采用 两 种 互信 
尽 的 方法 将 SAR ARS SPOT 图 像 进 行 配 准 . 文献 1 了] 采用 方位 互信 息 改进 了 原 有 的 互信 息 方 
法 ， 但 忽略 了 四 个 方向 上 的 空间 信息 . 文献 [8] 改进 了 方位 互信 息 度量 ， 用 开 窗 的 方法 计算 对 
比 度 图 像 ， 窗 口 包含 局 部 空间 信息 较 方 位 互信 息 方 法 有 效 ， 但 计算 量 大 . 

基于 特征 的 方法 并 不 直接 对 灰 度 级 图 像 进行 操作 ， 而 是 通过 提取 参考 图 像 和 待 配 准 图 像 上 
的 一 些 共 同 的 特征 ， 建 立 特 征 之 间 的 对 应 关系 ， 进 而 实现 图 像 的 配 准 . 因此 ， 特 征 提取 和 特征 
匹配 就 成 为 基于 特征 的 图 像 配 准 技术 中 的 两 个 关键 性 步骤 四 .常用 的 特征 有 点 特征 和 边缘 或 轮 
廓 特征 . 文献 [10,11] 中 基于 Canny 算 子 得 到 好 的 边缘 检测 结果 ， 但 SAR 图 像 存在 纹理 特性 很 
多 连续 轮廓 会 出 现 断 裂 现 象 ，LoG 算 子 则 会 产生 很 多 虚假 轮廓 .Lin2,13] 用 简单 的 轮廓 和 弹性 
轮廓 模型 分 别 对 光学 和 遥感 图 像 进行 配 准 . 简单 的 轮廓 算法 采用 链 码 相关 来 匹配 两 幅 图 像 的 轮 
廓 。 链 码 相关 方法 在 对 于 没有 明显 轮廓 的 图 像 配 准 问题 很 难 奏 效 ， 且 它 对 噪声 敏感 . 弹性 轮廓 
用 主动 轮廓 模型 ， 由 于 SAR 图 像 斑 点 噪声 干扰 ， 仅 用 轮廓 线 来 做 配 准 信息 量 少 ， 只 有 在 图 像 
差异 不 大 的 情况 下 才能 得 到 较 好 的 结果 . 
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针对 以 上 方法 用 于 SAR 图 像 配 准 的 不 足 之 处 ， 本 文 提 出 一 种 基于 主动 轮廓 分 割 区 域 
的 SAR 图 像 配 准 方法 . 该 方法 首先 利用 主动 轮 廊 模型 对 预 处 理 后 的 SAR. 图 像 进 行 分 制 ， 得 到 
分 割 后 的 若干 闭合 区 域 ， 用 区 域 特征 点 作为 特征 并 利用 区 域 匹 配 度 函数 作 度 量 求 得 粗 匹 配 下 的 
几何 变换 参数 ， 最 后 利用 互信 息 方法 在 最 佳 区 域 吃 配对 上 优化 参数 ， 并 将 图 像 进行 配 准 . 该 方 
法 克服 了 传统 边缘 检测 算 子 (如 Canny、Sobel 和 LoG 等 ) 受 SAR 图 像 斑点 噪声 干扰 较 大 的 问 
题 . 选取 区 域 对 作为 SAR 图 像 配 准 的 优势 在 于 : 其 一 ， 区 域 包 含 的 信息 较 轮 廓 的 多 ， 因 而 不 
会 因为 检测 轮廓 的 不 准确 导致 比较 差 的 配 准 结果 ; 其 二 ， 在 区 域 对 上 用 互信 息 估 计 参 数 计算 复 
杂 度 大 大 降低 ， 对 斑点 噪声 不 敏感 ， 用 于 SAR 图 像 配 准 效果 好 . 

本 文 后 续 内 容 安 排 如 下 : 第 2 节 给 出 算法 的 数学 模型 ， 第 3 节 介 绍 了 基于 分 割 区 域 
的 SAR. 图 像 配 准 的 计算 步骤 . 第 4 节 给 出 了 实验 结果 ; 最 后 一 部 分 为 结论 . 


2 主动 轮廓 与 互信 息 的 数学 模型 


2.1 主动 轮廓 模型 

主动 轮廓 是 由 Kassl4 等 人 于 1987 年 提出 的 . 其 基本 思想 是 以 构成 一 定形 状 的 一 些 控制 点 
为 模板 (轮廓 线 )， 通 过 模板 自身 的 弹性 形变 ， 与 图 像 局 部 特征 相 匹 配 达到 某 种 意义 下 的 协调 ， 
即使 得 菜 种 能 量 函 数 极 小 化 ， 完 成 对 图 像 分 割 ， 优 点 在 于 图 像 中 图 形 形 状 不 会 因 控 制 点 的 移动 
而 受到 影响 . 

用 一 组 向 量 序列 v(s) = (z(s),y(s)) 代表 一 个 轮廓 H3 ， 其 中 z(s) 和 w(s) 为 像素 的 坐标 ，1 < 
s< N，NN 为 像素 点 的 个 数 . 记 IT(zx,y) 为 图 像 函数 ， 能 量 函 数 定义 为 


E =X (a(s)Bi+6(s)E2 + 7(s)E3), (1) 


和 这 二 


1 N 
EE lew 
i= 
Eo = lus—1 一 2v; + vst+1|2, Es = [VI(z, y)|?, 


其 中 (1) 式 中 Qa(s)，B(s) 和 ~y(s) 为 权重 因子 . Bi 和 Es 代表 了 图 像 内 部 的 能 量 ，Ei 为 图 像 边 缘 
点 的 一 致 分 布 而 Bo 为 光滑 准则 函数 ;Es 为 图 像 的 梯度 . 给 一 个 初始 轮廓 ， 检 测 轮 廊 线 上 每 个 
点 相 邻 的 点 ， 当 所 检测 点 使 得 能 量 王 的 值 减 小 时 ， 相 应 轮廓 线 上 点 的 位 置 移 至 该 检测 点 上 . 和 连 
代 直 到 能 量 函 数 达到 最 小 ， 这 样 便 得 到 了 最 终 的 轮 廊 线 . 

2.2 ” 仿 射 变换 模型 及 区 域 匹配 度 函 数 

本 文选 用 仿 射 变换 模型 作为 参数 变换 模型 . 仿 射 变换 模型 


(x)-(" 2) ): (Ye) e 


其 中 (X,Y 了) 是 参考 图 像 中 4A 点 的 坐标 对 应 待 配 准 图 像 坐标 A 为 (Xo,Yo)， 尺 度 参 数 k = 
Vu? 十 v2， 旋 转角 度 为 


V T NT 
0 = arcta (=), 一 二 <0< 一 
ATRN 2 <p 


变换 模型 (2) 可 表示 为 4 = T(A’). 下 面 介绍 区 域 匹 配 度 函 数 C. 
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Bri, ro 分 别 是 参考 图 像 和 符 配 准 赂 像 中 的 分 割 区 域 ， 选 取 两 对 特征 点 对 带 入 (2) 式 即 可 
求 得 变换 水 数 工 (.)， 将 变换 函数 T(.) 作用 于 ro 得 到 rs = T(rz)， 分 别 取 矩形 Mo x No 和 x 
Ni， 将 ri 和 ra 包含 在 内 ， 令 


M = max( Mo, Mı), N = max(No, Ni). 


EXU FAAR 
1, (x, y) = Ty, 
fo(z,y) = (3) 
a | 0, (xy) gr 
1, (x,y) ars, 
;Yj 二 4 


Wo) FX ri 和 ra 的 匹配 度 函 数 为 


C(r1, r2) = 


上 式 第 二 个 等 号 由 (3), (4) BE, OP RR Mrs 重合 区 域 的 点 的 个 数 ， 分 母 根 式 中 为 区 
Sr) 和 rs 点 个 数 的 乘积 . 显然 C(ri,72) 的 范围 在 [0,1] 之 间 . 当 m 和 rm 是 一 个 匹配 区 域 ， 
且 满 足 = T(r2) 时 ， 函 数 C(r1,72) 达 到 最 大 值 1. 但 是 在 实际 匹配 中 ， 并 不 能 得 到 理论 
值 ]， 这 主要 由 于 两 幅 图 像 对 应 分 割 区域 本 身 可 能 并 不 是 完全 重合 的 . ARR r 和 rz 是 
匹配 区 域 ， 即 ri 和 rs 重合 部 分 较 多 ， 则 求 得 的 C(ri,72) 值 越 大 ， 为 一 方面 C(r1,72) 的 值 越 
大 ， 即 ri1 Arg BBM RS, WA Kr 和 7 是 匹配 区 域 的 概率 越 大 . 因此 本 文选 取 使 
得 C(r1,72) 达到 最 大 的 那 对 匹配 区 域 作 为 最 佳 区 域 对 作为 细 匹 配 ， 

2.3 互信 息 

互信 息 (Mutual information，MJ)5515 Æ Shannon 在 1948 年 关于 信息 论 的 论文 中 提出 来 的 
一 个 重要 概念 . 互信 息 是 用 来 度量 两 个 随机 变量 之 间 的 相似 程度 . 对 于 两 个 随机 变量 4, B, 
ETT BY 5 Sp E H(A), HH(B) 和 万 (4, B)， 概 率 分 布 分 别 为 Pa(a) 和 Pp(b)， 联 合 概 率 分 布 
为 Pa,p(a,b)， 则 其 互信 息 MI(A, B) 为 


Pap (a, b) 


Pa(a)Pp(d) 6) 


MI(A, B) = H(A) + H(B)H(A, B) = 》 Pas(a, b) log 


互信 息 是 一 种 图 像 相 似 性 度量 .在 基于 互信 息 的 图 像 配 准 中 ， 将 像素 值 看 作 是 某 个 连续 随机 变 
量 的 采样 ， 利 用 这 些 离散 采样 值 估计 出 单个 图 像 的 概率 密度 和 两 幅 图 像 的 联合 概率 密度 ， 再 求 
出 图 像 的 互信 息 . 两 幅 图 像 配 准 得 越 好 ， 则 它们 的 相关 性 越 大 ， 得 到 的 互信 息 值 也 就 越 大 ， 反 
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之 ， 两 幅 图 像 互 信息 值 越 大 ， 说 明 它们 之 间 的 相关 性 越 大 ， 则 图 像 配 准 得 越 好 .由 此 可 见 ， 基 
于 互信 息 量 测度 下 的 图 像 配 准 可 以 表示 为 


T*(-) = arg max MI(A,T(B)), (7) 


EH T*(-) ARRAS RM. 
互信 息 配 准 方法 目前 被 公认 为 是 配 准 精 度 和 和 鲁 棒 性 最 好 的 配 准 方法 之 一 . 


3 ”基于 分 割 区 域 SAR. 图 像 配 准 的 计算 步骤 


本 节 给 出 基于 分 割 区 域 SAR 图 像 配 淮 的 计算 步骤 : 

步骤 1 将 参考 图 像 与 待 配 准 图 像 分 别 进行 滤波 预 处 理 ; 

步骤 2 ”利用 主动 轮廓 模型 分 割 滤 波 后 图 像 ， 得 到 两 幅 图 像 的 分 割 区 域 ; 

PRS 提取 图 像 分 割 区 域 的 特征 ; 

步骤 4 利用 区 域 匹配 度 函 数 将 图 像 进行 粗 死 配 ， 得 到 匹配 区 域 及 初步 变换 参数 ; 

PRS 利用 互信 息 模 型 优化 参数 ， 得 到 细 匹 配 变 换 参 数 将 图 像 进行 配 准 . 

步骤 1 选用 高 斯 滤波 器 对 SAR 图 像 进行 滤波 ， 步 又 2 应 用 2.1 节 介绍 的 主动 轮 廊 模 型 将 图 像 
分 割 并 得 到 闭合 区 域 . 步骤 3 为 提取 分 割 区 域内 部 特征 过 程 ， 用 区 域 的 质心 以 及 Harris 角 点 检 
测算 法 提取 分 割 区 域 的 角 点 作为 特征 点 45j， 利 用 2.2 节 算法 求 取 粗 死 配 下 的 变换 参数 及 最 佳 
区 域 匹 配对 ， 得 到 最 佳 匹配 对 及 变换 参数 ， 再 利用 互信 息 方法 在 最 佳 匹配 区 域 上 优化 粗 配 准 
BR. 由 2.3 节 知 使 互信 息 达 到 最 大 求 得 的 配 准 参数 即 是 最 终 优化 参数 . 用 灰 度 直方 图 作为 图 
像 的 概率 密度 估计 ， 通 过 计算 两 幅 图 像 的 灰 度 联合 直方 图 来 估计 联合 概率 密度 . 优化 方法 采 
用 Powell 多 参数 优化 算法 ， 由 于 只 对 最 佳 区 域 对 实现 优化 算法 ， 这 样 有 效 地 减少 了 计算 量 . 


4 ”实验 结果 


4.1 SPOT MRF Landsat 图 像 配 准 结果 

图 像 来 自 文献 [13]， 图 1 是 利用 主动 轮廓 模型 分 割 两 幅 图 像 的 结果 . 图 1 中 (a) 和 (b) 若干 小 
区 域 未 检测 为 分 割 区 域 ， 这 是 由 于 实现 分 割 算法 时 剔除 了 小 区 域 ， 选 取 区 域 像素 个 数 大 于 图 
像 本 身 像素 个 数 1/100 (可 调节 ) 的 区 域 作为 分 割 区 域 . 图 2 给 出 了 分 割 区 域 特 征 与 配 准 结果 . 
选取 SPOT 图像 为 参考 图 像 ，Landsat 图 像 作为 待 配 准 图 像 进行 配 准 . 从 图 2 可 以 看 出 仅 用 
像 的 质心 作为 特征 点 匹配 图 像 的 不 足 之 处 在 于 对 应 的 分 割 区 域 并 不 一 定 完全 重合 ， 如 图 1 中 区 
Jat (c) 只 是 待 配 准 图 像 区 域 i) 中 的 一 部 分 ， 它 们 的 质心 在 几何 上 是 没有 对 应 关系 的 . 因此 本 文 
选取 角 点 特征 与 区 域 质 心 一 并 作为 粗 匹配 的 特征 更 加 合理 . 最 终 得 到 配 准 的 变换 参数 为 尺度 参 
数 为 

k = 1.371, @=17.9221°, AX =—59.81614, AY = 16.535. 


4.2 ”SAR 图 像 的 配 准 结果 

图 3 是 对 两 幅 SAR 图 像 进行 配 准 的 结果 ， 其 中 (b) 和 (e) 分 别 为 分 割 区 域 及 特征 ，(b) 图 为 
五 个 区 域 ， 而 (ce) 图 则 有 7 个 分 割 区 域 ， 但 两 幅 图 右 侧 的 四 块 对 应 区 域 都 被 检测 出 来 ， 这 样 对 
配 准 结果 并 无 影响 ， 因 为 只 需 检 测 出 一 对 匹配 的 区 域 即 可 完成 算法 .由 第 三 节 计算 步骤 求 得 图 
像 细 匹配 的 变换 参数 为 


k = 0.7234, 0 = 12.5077°, AX = 34.8267, AY = 88.5206. 
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图 4 是 对 男 两 幅 SAR 图 像 进 行 配 准 的 结果 ， 其 中 参考 图 像 (a) 大 小 为 100 x 100， 待 配 准 
图 像 (b) 的 大 小 为 128 x 128， 利 用 边缘 算 子 在 检测 这 两 幅 SAR 图 像 时 都 无 法 得 到 准确 的 轮 
廓 (c) 和 (d)，(e) 和 (f) 给 出 的 分 割 结果 和 特征 反应 了 本 文 算法 的 有 效 性 ， 最 终 得 到 图 像 细 匹配 
的 变换 参数 为 


k = 0.7813, @=13.1780°, AX =12.3850, AY = —9.2483. 


5 结论 


在 基于 特征 的 配 准 图 像 中 ， 轮 廊 方 法 的 优点 是 能 够 较 好 地 反应 图 像 本 身 的 信息 ， 计 算 复 杂 
度 低 . 但 是 传统 的 算 子 提取 SAR 图 像 边 缘 时 ， 提 出 的 轮廓 线 往 往 是 断裂 的 ， 甚 至 会 得 到 虚假 
轮廓 ， 直 接 用 轮廓 做 配 准 很 难 达到 满意 的 结果 . 本 文 利 用 主动 轮廓 分 割 图 像 ， 在 得 到 分 割 结果 
的 情况 下 ， 不 直接 提取 边缘 作 特 征 来 配 准 ， 而 是 把 得 到 的 区 域 作 为 特征 ， 先 进行 粗 匹 配 得 到 匹 
配 的 区 域 对 ， 再 利用 互信 息 得 到 细 配 准 的 结果 . 两 幅 SAR 图 像 配 准 和 遥感 图 像 配 准 的 结果 均 
表明 算法 是 有 效 的 .由 于 算法 本 身 是 基于 先 分 割 ， 再 配 准 ， 所 以 寻求 较 好 的 分 割 模型 作为 配 准 
前 提 是 笔者 今后 致力 研究 的 问题 . 
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图 1: (a): SPOT 图 像 ; (b): Landsat 图 像 ; (c)~(g): SPOT 图 像 的 分 割 区 域 ; 
(h)~(p): Landsat 图 像 的 分 割 区 域 
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图 2: (a): SPOT 图 像 分 割 区 域 特征 ; (b): Landsat 图 像 分 割 区 域 特征 ; 
(c): 变换 参数 作用 于 Landsat 图 像 的 结果 ; (d): 配 准 结果 





(d) (e) (f) 


图 3: (a),(d) 分 别 为 原始 SAR. 参 考 图 像 和 待 配 准 图 像 ; (b): SAR 参考 图 像 分 割 区 域 特征 ; 
(c): 变换 参数 作用 于 待 配 准 图 像 (d) 的 结果 ; (e): SAR 待 配 准 图 像 分 割 区 域 特征 ; 
(f): 配 准 结果 
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(e) (£) (g) (h) 


图 4: (a), (b) 分 别 为 原始 SAR. 参 考 图 像 和 待 配 准 图 像 ; (c): Canny 算 子 检测 待 配 准 图 像 (b) 的 结 


果 ; (d): LoG 算 子 检测 图 像 (b) 的 结果 ; (e): SAR 参考 图 像 分 割 区 域 特征 ; (f): SAR HAL 
准 图 像 分 割 区 域 特征 ; (g): 变换 参数 作用 于 待 配 准 图 像 (b) 的 结果 ; (h): 配 准 结果 
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Abstract: Due to the speckle imaging and the impact of the different conditions, the synthetic aper- 
ture radar (SAR) images may differ considerably under the same scene. However, using the method of 
the gray or edge feature, the results of SAR image registration may not be ameliorated. Therefore, we 
propose a method based on the segmentation-derived region to register the SAR image. After denoising 
the SAR images, the segmentation regions can firstly be obtained by using the active contour image 
segmentation. Moreover, we extract the feature of segmentation-derived region, using the function of 
regional the matching degree to carry out its rough match. Finally, by using mutual information, we 
optimize transformation parameters such that a coarse to fine image automatic registration can be 
implemented. The experimental results indicate that the algorithm is effective. 

Keywords: synthetic aperture radar; active contour; segmentation-derived region; mutual informa- 


tion; automatic registration 
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